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With a sample of 232 xlO 6 T(45) -> BB events collected with the BABAR detector, we study 
the decays of B mesons to pph final states, where h = iv + , K$, K*° or K* + . We report evi- 
dence for the B° — ► ppK*° decay, with a branching fraction (1.5±0.5(stat)±0.4(syst)) x 10 -6 , and 
for the B + — ► n c K* + decay, with the branching fraction of B(B + — > n c K* + ) x B(rj c — > pp) = 
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(1.57lQ'45(stat)lQ'3g (syst)) x 10 -6 , and provide improved measurements of the branching fractions 
of the other modes of this type. We also report the measurements of the charge asymmetry consis- 
tent with zero in the B + — > ppTx + , B° —> ppK*° and B + — > ppK* + modes. No evidence is found 
for the pentaquark candidate 6> + in the mass range 1.52 to 1.55 GeV/c 2 , decaying into pK° s , or the 
glueball candidate /j(2220) in the mass range 2.2 < m P p < 2.4GeV/c 2 , and branching fraction 
limits are established for both at the 10~ 7 level. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

This paper reports measurements of the branching 
fractions of the baryonic three-body decays B — » pph 
where h = ir + , K®, K*°, or K* + and their resonant sub- 
structures. This work is a continuation of our study of 
the B + — > ppK + decay reported in Ref . [1] . Observations 
of these decays were reported earlier in Ref. [2], except 
for the channel B° — > ppK*°, for which only an upper 
limit was obtained. 

These channels are interesting for the dynamical infor- 
mation in the distribution of the three final-state parti- 
cles and for the possible presence of exotic intermediate 
states, such as a pentaquark candidate 6 + (1540) in the 
m pK o spectrum [3] or an /j(2220) glueball candidate in 
the mpp spectrum [4]. 

The quark diagrams of the three-body baryonic B de- 
cays arc described in detail in Ref. [5]; only the dom- 
inant contributions are described below. In particular, 
the B + — > ppn + decay proceeds mainly through external 
and internal W— emission "tree" processes, while in the 
decay B° — > ppK ^ a virtual loop "penguin" process 
b — > sg is expected to be dominant. The B + — > ppK + ^ 
decay receives contributions from both the penguin and 
the doubly Cabibbo-Kobayashi-Maskawa-suppressed ex- 
ternal W— emission tree processes. 

An important feature of B — > pph decays is an en- 
hancement at low pp masses reported in Ref. [2], sim- 
ilar to that seen in several other baryonic decays of B 
mesons [6] and J/ip [7]. This enhancement might reflect 
short-range correlations between p and p in a fragmenta- 
tion chain [5, 8]. Alternatively, it could be a feature of a 
quasi-two-body decay in which the pp system is produced 
through an intermediate resonance [4], in particular the 
X (1835) state recently observed by BES [9]. Rosner [10] 
proposed distinguishing between the two mechanisms by 
studying the distribution of events in the Dalitz plot. As 
in the case of the resonance, p and p are produced inde- 
pendent from a hadron, then the distributions m ph and 
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mph should be identical. 

Ten experimental groups reported narrow enhance- 
ments near 1.54GeV/c in the invariant-mass spectra for 
nK + and pK° [11]. The minimal quark content of a state 
that decays strongly to nK + is dduus: therefore, these 
mass peaks were interpreted as a possible pentaquark 
state, called 0(154O) + . On the other hand, a number of 
experiments that observe strange baryons with mass sim- 
ilar to that of the 6(1540)+ [e.g., 71(1520) -» pK~ ] see 
no evidence for the O(1540) + [12, 13]. As suggested in 
Ref. [3] we search for a pentaquark baryon candidate + , 
in the m pK o mass distribution of B° — > ppK° decays. 

A narrow state /j(2220) with a width of 23MeV was 
reported in the KK spectrum by the MARK III ex- 
periment [14] and later seen in K + K~, K°K , 7r + 7r~ 
and pp by the BES [15] experiment. However, its non- 
observation in quite a few pp annihilation modes, in par- 
ticular in pp -> /j(2220) KK [16], has led to doubts 
about its existence [17]. One theoretical conjecture [4] 
suggests a possible presence of the /j(2220) resonance in 
baryonic B decays. 

Direct CP violation is observable as an asymmetry in 
yields between a decay and its CP conjugate when at 
least two contributing amplitudes carry different weak 
and strong phases. Following the observation of the di- 
rect CP violation in B° — > K + ir~ [18] transitions and 
its non-observation in B + — > K + ir° [19], it is interesting 
to study the charge asymmetry in the B — > pph system. 
Three-body baryonic B decays to pph occur through two 
different mechanisms (penguin and tree), which carry 
different weak phases and, in general, different strong 
phases. Although negligible direct CP violation asymme- 
try is measured in the B + — > ppK + mode [1, 2], recent 
theoretical calculations predict a significant asymmetry 
in the B+ -> ppK*+ mode [20]. 

We use events with the same final state particles to 
isolate the decay B° — > A+p, A+ — > pK ^ and decays 
B — > X C ch, for X C c = rj c and J/ip, with a subsequent 
decay to pp. Charge-conjugate reactions are included 
implicitly throughout the paper. 

In Sec. II, we describe the BaBAR detector. In Sees. 
Ill, IV, and V we discuss the reconstruction of the 
B-meson candidates, the branching fraction (B.F.) ex- 
traction and the systematic uncertainties, respectively. 
In Sec. VI we present our results on the B — ► pph, 
B — ► X c -h and B° — ► A+p branching fraction measure- 
ments, charge-asymmetry measurements, and searches 
for O(1540) + pentaquark and /j(2220) glueball states. 
We discuss the B — > pph decay dynamics and summarize 
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in Sees. VII and VIII, respectively. 



II. BABAR DETECTOR AND DATA SAMPLE 

This measurement is performed with a data sample of 
232 xlO 6 T(4S) -> BB decays corresponding to an inte- 
grated luminosity of 210 fb _1 , collected at the T(4S) res- 
onance ( "on- resonance" ) with the BABAR detector at the 
PEP-II e+e - storage ring. An additional 21 fb -1 of data 
( "off- resonance" ) , collected about 40 MeV below the reso- 
nance, is used to study the backgrounds from light-quark 
and cc production. The BABAR detector is described in 
detail elsewhere [21]. Charged-particle trajectories are 
measured by a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH), both operating in a 1.5- 
T solenoidal magnetic field. A Cherenkov radiation de- 
tector (DIRC) is used for charged-particle identification. 
The CsI(Tl) electromagnetic calorimeter detects photon 
and electron showers. The BABAR detector Monte Carlo 
simulation based on GEANT4 [22] is used to optimize 
selection criteria and to determine selection efficiencies. 



III. B MESON RECONSTRUCTION 

The B + and B° decays are reconstructed in the fol- 
lowing modes: B + — > ppir + , B° —> ppK* a , B° —> ppK® 
and B + — > ppK* + . Charged tracks not coming from K® 
must originate from the interaction point, have trans- 
verse momentum greater than O.lGeV/c, and at least 
twelve DCH hits. Charged pions, kaons and protons 
are identified by the average energy loss (dE/dx) in the 
tracking devices and by the pattern of Cherenkov pho- 
tons in the DIRC. The — ^ 7r + 7r~ candidates are re- 
quired to have a two-pion invariant mass within 8 MeV/c 2 
of the nominal K® mass [12] and a cosine of the angle 
between the line connecting the B and decay ver- 
tices and the momentum greater than 0.999 in the 
laboratory frame. K* + candidates are reconstructed in 
the K* + — > Kgn + decay channel, and K*° candidates 
are reconstructed in the K*° — ► K + n~ mode. Candi- 
dates whose mass is within 80 MeV/c 2 of the nominal K* 
mass [12] are selected as K* candidates, and the ones 
with mass 160 — 240MeV/c 2 away from the K* mass are 
chosen as the K* sideband sample. The three daughters 
in the B decay are fit constraining their paths to a com- 
mon vertex. To suppress background, we impose a cut 
requiring a probability greater than 10~ 4 on the each of 
the K° and the B candidates' vertices. 

Two kinematic constraints of B-meson pair-production 
at the T(45) are used to characterize the B candidates: 
the beam-energy-substituted mass ities = [(E^ m /2 + 
Po ' Pb) 2 /Eq — P^] 1 / 2 and the energy difference AE = 
E* B — E cm /2, where E cm is the total center-of-mass en- 
ergy, Eg is the center-of-mass energy of the B me- 
son candidate (the asterisk denotes the center-of-mass 
frame throughout the paper), and (Eq,p ) is the four- 



momentum of the initial state and p B is the B momen- 
tum in the laboratory frame. For signal B decays the 
AE distribution peaks near zero with a mode-dependent 
resolution of 11 — 17 MeV, while the tues distribution 
peaks near the B meson mass with a mode-dependent 
resolution of 2.5 — 3.0MeV/c 2 , as determined in Monte 
Carlo simulation. 

We select events with m ES > 5.22 GeV/c 2 (5.25 GeV/c 2 
for the B+ -> ppn+ mode) and \AE\ < 0.10 GeV 
(0.15 GcV for the B Q -> ppK° mode). We make these 
selections quite loose since these two variables are used 
in a maximum likelihood fit to extract the signal yield. 
Of the candidates passing that loose selection, only one 
candidate is chosen for each event, selecting the one with 
the highest B vertex probability, or, if a K® is present 
in the decay chain, the highest value of the product of 
the B and the K® vertex probabilities. To improve the 
resulting mass resolutions, after selecting the B candi- 
dates we perform a kinematic fit fixing the mass of each 
B candidate to its known value [12] and its energy to a 
half of the total center-of-mass energy. 

The background is dominated by random combinations 
of tracks created in e+e - — ► qq (where q is u, d, s,c 
quarks) continuum events. These events are collimated 
along the original quark directions and can be distin- 
guished from more spherical BB events. We construct a 
Fisher discriminant (J 7 ) [23] as a linear combination of 
the following four event-shape variables: 

1. cos9^ hr , the cosine of the angle between the thrust 
axis of the reconstructed B and the beam axis; 

2. cos#* lom , the cosine of the angle between the mo- 
mentum of the reconstructed B and the beam axis; 

3. and 4. the zeroth and the second Legendre poly- 
nomial momentum moments, Lq and Li. They 
are defined as follows: Lq = J2i \Pi \ an d L 2 = 
Y\ |p*|[(3cos 2 e* hrB i -l)/2], where p* are the mo- 
menta of the tracks and neutral clusters not associ- 
ated with the B candidate and 0^ hrB . is the angle 
between p* and the thrust axis of the B candidate. 

The coefficients for the Fisher discriminant variables are 
determined for each of the modes separately by max- 
imizing the separation between the means of the sig- 
nal and background Fisher discriminant distributions ob- 
tained from the signal Monte Carlo samples and the off- 
resonance data events, respectively. 

There is also a background for the three-body decays 
B — > pph from B mesons decaying into the same final 
states as the signal, such as B — > X c -h decays, where 
X C c = rjc, J/ip, ip' ',XcO,Xci,Xc2 decaying to pp and the 
B° — ► A+p decay. This background, which is of interest 
in its own right, is discussed in the next section. Other 
backgrounds from B decays are negligible. 
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IV. BRANCHING FRACTION EXTRACTION 

We perform an unbinncd extended maximum likeli- 
hood (ML) fit to extract the signal yields. The variables 
toes, AE and T are used as discriminating variables to 
separate signal from background. 

To estimate the contribution from B decays that pro- 
ceed through b — > c transitions leading to the pph final 
state, such as rj c , J/ip — > pp and A£ — > pKg/K* a , the 
maximum likelihood fit is performed in three distinct re- 
gions: 

1. The main "charmonium" region which includes the 
■q c and J/ip resonances: 2.85 <m pp < 3.15GeV/c 2 
(for h = K°/K*°/K* + ). 

2. The "charm" (A£) region delimited by \m p h — 
2.3 GcV/c 2 | < 0.1 GeV/c 2 where h = K°/K*°, and 
\m pp — 3 GeV/c 2 | > 0.15GeV/c 2 (this prevents over- 
lap with the charmonium region). Note that for 
h = Kg both m p h and m p h combinations are con- 
sidered. 

3. The "all-other" region. Significant background 
from slow pions in the B + — ► ppn + mode leads to a 
slight difference of the Fisher discriminant shape at 
low and high m pp . To reduce the sensitivity to that 
correlation we perform the fit for this mode in two 
regions: m pp < 3.6GeV/c 2 and m pp > 3.6GeV/c 2 . 

In the charmonium region, in addition to the three 
variables described above, the m pp variable is used to 
distinguish between the non-resonant signal and contri- 
butions from to, and J/ip. In the charm region the corre- 
sponding additional variable is m p h- 

The data sample is assumed to consist of two compo- 
nents: signal events, including B meson decays to X c -h 
and A£p, which have the same final state particles as the 
signal, and combinatorial background events. The latter 
are due to random combinations of tracks from both con- 
tinuum and BB events. For the fit in the charmonium 
region, rj c and J/ip signal-like components are included 
in the fit, a A£ component corresponding to the overlap 
of charm and charmonium regions is not included here 
but its contribution is taken into account in Sec. VIA. 
For the charm region, a component is included in 
the fit as both signal and background contributions. For 
the decay B + — > ppn + , an additional B + — > ppK + com- 
ponent is considered, to account for the B + — ► ppK + 
events present in the fit region because of kaon/pion mis- 
identification. The signal component is split into cor- 
rectly reconstructed true events and mis-reconstructed 
events, so-called self-cross- feed (SCF). The SCF events 
are signal events in which one or more of the tracks from 
the signal side are lost and a track from the other B 
decay is included in the reconstruction. The fraction of 
SCF events is determined with a B — > pph Monte Carlo 
sample and varies from 0.5% for the B° — ► ppK® mode 
to 5.6% for the B+ ppK*+ mode. In the B° -> ppK° s 



mode, the SCF signal events are indistinguishable from 
the combinatorial background and no separate SCF sig- 
nal component is used in the maximum likelihood fit. 

In the maximum likelihood fit, each component is 
modeled by the product of probability density functions 
(PDF) of the following variables, which are assumed to 
be uncorrelated for all components: 

• Toes, AE, T and m pp in the charmonium region: 

n = V x (m ES ,AE,F,m pp ) (1) 

• toes, A.E and T and m p h in the charm region: 

VI = V x {m^,AE,T,m ph ) (2) 

• toes, AE and T in the all-other region: 

V%=V x (m ES ,AE,F) (3) 

where x indicates the event source. In the all-other re- 
gion there are two components: x is either signal (s) or 
background (b) , in the charmonium region there are four 
components x = s, i) c , J/ip and b; and in the charm re- 
gion there are three components x = s, A+ and b. 
The likelihood is given by 

i— 1 x 

where r corresponds to a fit region, N r is the total num- 
ber of events in that region and ls * ne sum over au 

X 

the corresponding fit components in a given region. Then 
the total number of events in the charmonium region is 
N\ = Nf + + N J ^ + iVf, in the charm region is 
N 2 = N$ + N A i + AT|, and in the all-other region is 
N 3 = N§ + V|, where N s , N^, N J ^, N A t and N b are 
the number of non-resonant signal, r) c , J/ip, signal 
and combinatorial background events, respectively. V x 
is the PDF for a corresponding component x, e.g. V^. 
is the background PDF for region r. The signal PDF 
V s r = P* rue + /scF^ r SCF consists of the PDFs of the 
true and SCF signal events, respectively, with the cor- 
responding fraction of SCF events /scf as determined 
from B — > pph signal Monte Carlo and fixed in the fit. 

The paramctrization of the B — ► pph true signal, SCF, 
r] c , J/ip, A£ and combinatorial background components 
is summarized in Table I. 

All parameters in the B — > pph true and SCF signal 
PDFs are obtained from corresponding B — > pph signal 
Monte Carlo simulations. The B — > pph signal Monte 
Carlo events are simulated according to a three-body 
phase-space decay expectation, which has been shown 
in Refs. [1, 2] not to reproduce the data. In order to im- 
prove the parametrization for the B — > pph signal PDFs 
the signal Monte Carlo events are re-weighted accord- 
ing to the m pp distributions from B + — > ppK + given in 
Refs. [1, 2]. 
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TABLE I: The PDF parametrization of B — ► pp?t signal, self-cross-feed (SCF), rj c , J/tp, A~£ and combinatorial background 
contributions. We use the following notation: "BifGauss" for a Bifurcated Gaussian, "G" for a Gaussian (2xG are two 
Gaussian distributions with common mean) , "Voigtian" for a convolution of a Breit-Wigner distribution and a sum of Gaussian 
distributions with common mean, "ARGUS" for a threshold function [24], "norm" for normalization and "ratio" is a ratio of 
the normalizations of a linear to a Gaussian contributions. The parameters floated for each of the functions are specified in the 
brackets. 



Component rriEs 


AE 


T 


Ttlpp 


m p h 


B u - ppK% 



signal BifGauss (mean) 

J/tp same as signal 

rj c same as signal 

A£ same as signal 

combinatorial ARGUS (slope) 



G+G (narrow mean) G+G constant (norm 

same as signal same as signal 2xG (narrow a 

same as signal same as signal 2 x Voigtian 

same as signal same as signal — 

linear (all) G+G (means) linear 



constant (norm) 



2xG 

linear+2xG (ratio) 



signal BifGauss (mean) 

J/tp same as signal 

rj c same as signal 

SCF ARGUS 
combinatorial ARGUS (slope) 



G+G (narrow mean) G+G constant (norm 

same as signal same as signal 2xG (narrow a 

same as signal same as signal 2 x Voigtian 

linear G+G constant 

linear (all) G (all) linear 



ppK w 



signal BifGauss (mean) 

J/rp same as signal 

rjc same as signal 

At same as signal 

SCF ARGUS+G 
combinatorial ARGUS (slope) 



G+G (narrow mean) G+G constant (norm 

same as signal same as signal 2xG (narrow a 

same as signal same as signal 2 x Voigtian 

same as signal same as signal — 

quadratic G+G constant 

linear (all) G (all) linear 



constant (norm) 



2xG 



linear+2xG (ratio) 



• pp-K^ 



signal BifGauss (mean) 

SCF ARGUS 
ppK + BifGauss 
ppK + SCF ARGUS 
combinatorial ARGUS (slope) 



G+G (narrow mean) G+G 
linear G+G 



G+G 
linear 
linear (all) 



G+G 
G+G 
G+G (all) 



The toes, A£, and T PDFs for the r} c , J/ip and A+ 
contributions are taken to be the same as the correspond- 
ing PDFs for the B — > pph signal. The known masses of 
the J I ip and A£ and the known mass and width of the 
rjc are used in the fits. All the other parameters that 
are not floated are fixed to the values obtained from the 
corresponding Monte Carlo simulation. 

All the parameters for the B + — > ppK + background 
events in the B + — ► ppir + mode are fixed in the fit to 
the values obtained from the B + — > ppK + Monte Carlo 
events re-weighted as described above as well as the ex- 
pected amount of the charmonium events with ppK + fi- 
nal states. 

For the combinatorial background events all the pa- 
rameters fixed in the fit are obtained from the on-T(45 f ) 
resonance data toes sideband (toes < 5.26GeV/c 2 ). 

The parametrization of the signal and combinatorial 
background PDFs does not vary with the fit region. 

The remaining floating parameters of the PDFs are de- 
termined by a maximum likelihood fit to the data. The 
fit uses the kinematic variables described above but is in- 
dependent of the location of the event in the Dalitz plot. 
This fit not only determines the various parameters of 



the PDFs, but also the number of signal and background 
events, and the covariance matrix for these event num- 
bers. 

Once the maximum likelihood fit provides the best 
estimate of the PDF parameters, we use a weighting 
technique [25] to measure the branching fraction and re- 
construct efficiency-corrected mass distributions. This 
method allows us to take into account the dependence of 
the efficiency on the position of a candidate in the Dalitz 
plot. Wc assume that the distributions in the variables 
toes, AE 1 , and T are uncorrelated with the location in 
the Dalitz plot. 

The event-dependent weight to be signal, Wl (also 
known as the s Weight in the sPlots method [25]), is com- 
puted from the PDFs: 

s ' U 

where V* J is the value of the PDF for the compo- 
nent x (x = signal, background, charmonium-signal, 
charmonium- r] c , etc.) in the event j for the fit region 
r (as defined in Equations 1-3), and V r (s,x) is the co- 
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variance between the number of signal events iV* and 
the number of events for the component x, , in the fit 
region r defined by: 



V r (s,xY 



d 2 (-lnC) 



(6) 



The sum of W| over all events j is just N s = Nf + iV| + 
A| , while the sum of W| over all events in a small area in 
phase space gives, in the limit of high statistics, the cor- 
rect distribution of the signal in that area while preserv- 
ing the total signal yield. These weights (W|) allow opti- 
mal discrimination between signal-like and background- 
like events [25] and serve as a tool to reconstruct the 
resulting signal distributions. 

The resulting branching fraction is then calculated as 
follows: 



N 



N BB ■ £ .i ■ #sub' 



(7) 



N B g is the total 



where the sum is over all events j 
number of BB pairs, B su b is the product of the branch- 
ing fractions in the sub-decays and Sj is the reconstruc- 
tion efficiency of the event j. The reconstruction effi- 
ciency depends on the position of the candidate on the 
Dalitz plane and is obtained from the corresponding sig- 
nal Monte Carlo simulations. The statistical error on the 
branching fraction is given by 



AB 



{wi) 2 

71 



(8) 



V. SYSTEMATIC ERRORS 

The contributions to the systematic errors on the mea- 
sured branching fractions and charge asymmetries are 
summarized in Tables II and III. 

The T(4S) is assumed to decay equally to B°B° and 
B + B~ mesons. Incomplete knowledge of the luminos- 
ity and cross-section leads to a 1.1% uncertainty in 
all branching fraction measurements. Chargcd-tracking, 
particle-identification (PID) and reconstruction stud- 
ies of the data lead to small corrections applied to each 
track in these simulations. Limitations of statistics and 
purity in these data-Monte Carlo comparisons lead to 
residual systematic uncertainties. Limitation of Monte 
Carlo statistics employed in each analysis contributes a 
small uncertainty. The effects of binning on the efficiency 
are estimated by studying the variation in the resulting 
signal yield due to changes in the chosen bin size. 

A large control sample of B^ J/ip{e + e~, n + [i~)h is 
studied separately in data and Monte Carlo simulations 
to understand uncertainties arising from the choice of ver- 
tcxing cuts and from the parametrization of the PDFs 



for AE, ities an d T (this uncertainty is labeled "Pre- 
selection" in Table II). 

The uncertainty due to possible correlation between 
the fit variables is estimated by performing fits to a 
number of Monte Carlo experiments that consist of fully 
simulated signal events embedded in parametrized back- 
ground simulations. This so-called "Fit Bias" uncer- 
tainty in Table II) is on the order of a few percent. 

Where the Monte Carlo values are used to fix signal 
shape parameters in a fit, we obtain the uncertainty on 
the signal yield due to a change in each of the PDF pa- 
rameters by doing a number of toy Monte Carlo exper- 
iments with parametrized signal and background distri- 
butions. In each a fit is performed with the original value 
of the PDF parameter and the value shifted by its uncer- 
tainty, which is obtained from the B — > pph signal Monte 
Carlo simulation. The resulting uncertainties are added 
together in quadrature, taking into account the correla- 
tion between the PDF parameters, to obtain the total 
error for each PDF. The procedure is repeated using the 
correlation matrices between the variables to obtain the 
total error for each fit component, which is then added in 
quadrature to give the total error on PDF parametriza- 
tion. 

In a similar fashion, different fit ranges for charmo- 
nium and charm regions are employed and the resulting 
variation of the fit yields is taken as a systematic error. 

The potential correlation of the fit variables with their 
location on the Dalitz plot would slightly violate the main 
principle of the sPlot method. To reduce the sensitivity 
to that correlation the fit is performed in one, two and 



for 



pp 



regions for all the modes. The resulting vari- 



ation of the branching fraction is taken as a systematic 
error. Note that we constrained the B + — > ppK* + fit 
to make the signal component non-negative (in the char- 
monium region), which results in additional systematic 
error on the fit region. All the uncertainties above are 
summed in quadrature to give the "Fit Region" error in 
Table II). 

Branching fraction uncertainties [12] on B(B + — > 
Xh) x B(X — > pp), where X = x c [o,i,2]> V*' an d h — 
K^, K*° or K* + , and on the sub-branching fractions of 
K° and K* , which affect the B — > pph branching fraction 
measurements, are given in the "B Bkg / B.F. errors" line 
in Table II. 

For modes that contain K* mesons, the non-resonant 
Kir background is obtained by performing a maximum 
likelihood fit and extracting the branching fraction in 
the K* sideband region. The non-resonant Kit con- 
tributions to the branching fraction values, based on 
the K* sideband data, for the the B + — ► ppK* + 
and B° -> ppK*° modes are (0.34±0.74) x 1(T 6 and 
(0.23±0.30)xlCT 6 , correspondingly. The X*-sideband 
signal yields for J/^jK* and J/xjjK*+ arc 11.5±|;| and 
1.2^09 events; for i] c K*° and rj c K* + arc ± 5.7 and 
1.8i| o events; for the B -> A+p in the B° -> ppK*° 
mode 4.5^2 7 events. As no significant non-resonant Ktt 



9 



TABLE II: Systematic errors (in percent) and efficiency corrections for all modes. The efficiency correction is a weight applied 
event-by-event to signal efficiency to account for residual differences between signal Monte Carlo and data. 



Error Source 


ppK 1 ' 


ppK* + 


ppK* u 


ppn + 








pAtipK" 


) pA+( P K* u ) 


e+ 


BB counting 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


PID efficiency 


1.2 


1.6 


3.4 


2.1 


1.2 


1.6 


2.9 


1.6 


3.6 


-5.8/+6.6 


Track reconstruction 


1.6 


2.4 


3.6 


2.4 


1.6 


2.4 


3.2 


1.6 


3.2 


1.6 


Kg reconstruction 


2.8 


2.9 


n/a 


n/a 


1.1 


1.2 


n/a 


2.1 


n/a 


1.6 


Monte Carlo Statistics 


2.0 


-3.4/+4.4 


2.5 


2.6 


0.3 


1.3 


0.8 


1.2 


1.0 


0.3 


Dalitz plot binning 


2.1 


3.5 


1.4 


2.0 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


Pre-selection 


0.5 


0.8 


0.4 


0.3 


0.5 


0.8 


0.4 


0.5 


0.4 


0.5 


Fit Bias 


5.0 


3.0 


3.6 


1.1 


3.4 


3.4 


1.5 


20.0 


11.0 


n/a 


PDF Parametrization 


3.0 


-3.7/+3.2 


-1.9/+2.2 


6.2 


1.1 


-3.3/+3.8 


3.0 


1.3 


0.6 


n/a 


Fit Region 


5.7 


17.0 


6.5 


13.0 


5.0 


1.8 


3.0 


5.0 


3.4 


n/a 


B Bkg / B.F. errors 


1.0 


0.8 


2.9 


0.8 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


Non-resonant Kn 


n/a 


15.4 


25.2 


n/a 


n/a 


-21.9/+22.7 


15.5 


n/a 


-42.7/+46.9 


n/a 


Total(%) 


9.5 


24.4 


27.1 


15.2 


6.7 


-22.8/+23.6 


16.8 


21.0 


-44.5/+48.6 


-6.3/+7.1 


Pre-selection correction 
Kg correction 


0.976 
0.981 


0.967 
0.980 


0.968 

n/a 


0.983 
n/a 


0.976 
0.972 


0.967 
0.966 


0.968 
n/a 


0.969 
0.976 


0.968 
n/a 


0.976 
n/a 



TABLE III: Systematic errors for the charge asymmetry measurements. 



Error 


ppK*+ 


ppK* 


Pp-K + 






J/tpK*+ 


J/^K* l> 


Track Reconstruction 


0.002 


0.004 


0.002 


0.002 


0.004 


0.002 


0.004 


PID Efficiency 


0.003 


0.015 


0.001 


0.001 


0.002 


0.002 


0.001 


Monte Carlo Statistics 


0.012 


0.057 


0.021 


0.017 


0.004 


0.002 


0.016 


Fitting 


0.050 


0.020 


0.040 


0.066 


0.040 


0.018 


0.007 


Total 


0.05 


0.06 


0.05 


0.07 


0.04 


0.02 


0.02 



background is seen, we do not perform the sideband sub- 
traction, but rather add these background contributions 
and their statistical uncertainties in quadrature to give 
the "Non-resonant Kn v systematic uncertainties listed in 
Table II. 

All the uncertainties described above are added in 
quadrature in the "Total" line in Table II. The final 
branching fraction value obtained from Eq. 7 is scaled by 
Pre-selection and 11 K^" corrections. 



VI. RESULTS 
A. B — > pph Branching Fraction Measurements 

The event yields from the maximum likelihood fit are 
presented in Table IV, while Fig. 1 shows projections of 
the corresponding three-dimensional PDFs on the AE 
and mEs axes. The measured branching fractions cal- 
culated from Eq. 7 (scaled by Pre-selection and K° cor- 
rections) are summarized in Table V, where the effective 
efficiency for the modes observed is calculated accounting 
for the correct event distribution in the Dalitz plot. 

Our fitting method removes all B — > r\ c h and B — > 
J/iph (except in the B + — > ppn + mode) contributions 
and most of the B° — > A+p contributions. There is 



still some remaining B° — > A+p background contribu- 
tion from the B — > A+p events in the charmonium re- 
gion. Knowing the relative efficiencies of B° — > A+p 
Monte Carlo events inside and outside the charmonium 
region allows us to calculate the remaining B a — ► A+p 
background contribution to be (0.04 ± 0.02) x 10~ 6 and 
(0.06 ± 0.03) x 10~ 6 for the B° -> ppK and the B Q -> 
ppK* modes, respectively. 

The remaining unknown background comes from the 
B — > Xcoh events. It is possible to estimate the 
B° — > XcoK branching fraction using the correspond- 
ing B + —> Xd)K + branching fraction measurement. Be- 
cause of isospin symmetry one would expect ratios of the 
charged and neutral B mesons decaying into Xco an d Xci 
to be equal. Thus we estimate [12]: B(B° — > XcoK ) « 
B(B+ - Xc0 iT +) gfeg^l = (i.e+o.5) x 10 - 4 x 

salo 7 = (°- 12 ± °- 04 ) x 10_3 - This number needs to 
be multiplied by the B(xc0 -> PP) = (0.22 ± 0.03) x 10~ 3 
and results in the expected contributions to the branch- 
ing fraction from this mode of (0.026±0.009) x 10~ 6 . The 
resulting contribution to the absolute systematic error on 
the B background is 0.009 x 10~ 6 . The contribution from 
the B — > Xcoh events is ignored for other modes [26]. 

From Rcf. [12] we estimate the contribution of B + — > 
J / xpir + to the B + — > ppir + mode and the remaining char- 
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TABLE IV: Summary of the resulting yields from the maxi- 
mum likelihood fit for all modes. 



Region 


charmonium 


charm All-Other 


Total 


Component 




B° -» ppR's mode. 




Signal 

Vc 


17± 9 8 
23+? 


3l* 7011? 


90t\t 
23 +8 7 


J/i> 
At 


53± 8 7 


6-8±|;| 


531? 


Background 


1152±34 


1096±33 14769±122 


17017±131 



B — > ppK* mode. 



Signal 


0±9 




JZ -10 


52±14 


Vc 


12.31 4 ;* 




12.31*:* 


J/i> 


34±6 






34±6 


Background 


766±28 




10063±101 


10829±105 






B° - 


* ppK* mode. 




Signal 

Vc 


O+10 

8 -8 

371^° 




4l 6 4 5011* 


63±il 
371^° 


j/tP 


106±11 






106±11 


At 




12.3 


+4.9 
-4.1 


12.31*;? 


Background 


2207±47 


1971±45 26312±163 


30490±176 






B+ - 


-> ppn + mode. 




Signal 






185±28 


185±28 


ppK + 






157±30 


157±30 


Background 






90438±305 


90438±305 



22 5.24 5.26 5.28 5.3 




0.05 0.1 
A E(GeV) 



25 5.26 5.27 5.28 5.29 5.3 
m ES (GeV/c 2 ) 



FIG. 1: Distributions (points with error bars) of AE and 
toes for data candidates passing loose selection in the B° — > 
ppK (a)-(b), B+ -> ppK*+ (c)-(d), 5° -> ppK* (e)-(f), 
B + — > ppir + (g)-(h) modes. The superimposed solid curves 
represent the projections on the three-dimensional-fit PDFs 
on the respective axis, and the dashed line is the fitted signal 
contribution. Note that the excess of events on the left of the 
signal peak in (g) corresponds to the B + — > ppK + final states, 
where K + is mis-identified as ir + . The following loose selec- 
tion is applied: (a) m E s > 5.275 GeV/c 2 and T > 0.2; (b) 
\AE\ < 0.02GeV/c 2 and T > 0.2; (c) m ES > 5.275 GeV/c 2 
and T > 0.2; (d) \AE\ < 0.02GeV/c 2 and T > 0.2; (e) 
m ES > 5.275 GeV/c 2 and T > 0.6; (f) \AE\ < 0.015 GeV/c" 



and T > 0.6; (g) m ES > 5.275 GeV/c 2 
\AE\ < 0.015 GeV/c 2 and T > 1.5. 



and T > 1.5; (h) 



monium contributions (B — + Xd/i and -B ip(2S)h) for 
all the other modes. Since only an upper limit exists 
for the branching fraction of the B — > XrJih mode, it is 
not subtracted but taken as a one sided contribution the 
to systematic uncertainty. The total expected B back- 
ground contribution is quoted in the line "B bkgr. B.F." 

The line "Final B.F." in Table V summarizes the values 
of the charmless and charmoniumlcss B — > pph branching 
fractions after B background subtraction. 

We report the first evidence for the B° — > ppK* de- 
cay with a significance of 3.2cr (including systematic un- 
certainties). The statistical significance a throughout 
the paper is taken as \J — 2 ln(£(0)/£ m ax), where £(0) 
is the likelihood of the fit assuming zero signal events 
and Umax is the likelihood obtained in the full fit. For 
the B° — > ppK* decay the C(Q) is taken to be not at zero 
signal events, but at the expected number of the B back- 
ground events as discussed above. To obtain the value 
of significance including the systematic uncertainties the 
likelihood function is smeared with a Gaussian distribu- 
tion which has a width of the corresponding systematic 
uncertainty. 

The measurements of branching fractions for the B — > 
pph modes from Ref. [2] and this work are summarized in 
Table VI and compared to those of the two-body mesonic 
modes from Refs. [12, 27, 28]. 

The branching fractions are approximately two times 



smaller for the B + — > ppK* + and the B^ 



ppn n 
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TABLE V: Summary of the non-charm, non-charmonium B —* pph branching fraction calculation including B background 
subtraction. Note that effective efficiency includes the sub-decay branching fractions as well as the correct event distribution 
in the Dalitz plot. 





ppK u 


ppK* + 


ppK* u 


ppiY + 


Sum of Signal Weights 


181±29 


52±14 


63±17 


185±28 


Effective Efficiency, % 


24.5 


4.1 


16.0 


44.6 


B.F. from Eq. 7 (10~ 6 ) 


3.17±0.53 


5.45±1.49 


1.70±0.45 


1.79±0.29 


B bkgr. B.F. (1(T 6 ) 


0.22±0.03 


0.17±0.04 


0.23±0.05 


0.10±0.01 


Final B.F. (10 _b ) 


3.0±0.5±0.3 


5.3±1.5±1.3 


1.5±0.5±0.4 


1.7±0.3±0.3 



modes, when compared to the Belle measurements [2], 
bringing the branching fraction of the B + — > ppK* + 
mode below that of the B + — > ppK + mode and more 
in line with theoretical predictions [5]. However, the two 
experiments are in agreement within their errors. 

Since the virtual loop "penguin" process b — ► sg pre- 
serves isospin we would naively expect the ratio of the 
rates for B + — > ppK + and B ppK° to be unity as it 
is in two-body mesonic modes B — > p°h and B —> ir°h, 
but it is closer to two (see Table VI). This could be ex- 
plained by absence of the external W— emission Feynman 
tree diagram for the neutral B mode. However, if this 
tree diagram were important, we would expect a much 
larger rate for B + — > ppn + than for B + — ► ppK + , in con- 
tradiction with the data. The B + — > ppK* + branching 
fraction is also larger (by a factor of three) than that of 
B° — > ppK*° similar to the pattern suggested by the data 
for decays to p°K* and tt°K*. The B — > ppK* modes are 
consistently smaller than the B — > ppK modes in both 
the charged and neutral cases. This seems to be the case 
for the B — > iz°h modes as well, but not for the B — ► p°h 
modes. The B + — > ppir + branching fraction is lower than 
that of the B + — > ppK + mode as expected because the 
b — > u transition at tree level is suppressed compared to 
the b — > s penguin. This is similar to what is observed in 
the B — > n h modes but contrary to what is observed in 
B -> p°h. 

Overall, the theoretical calculations of the baryonic B 
decays are not very certain and the current measurements 
of the branching fractions of all four B — > ppK modes are 
a further challenge to our understanding. 
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FIG. 2: The results of the maximum likelihood fit for ?bes 
> 5.27GeV/c 2 and \AE\ < 0.03 GeV for the B+ -> ppK* + 
mode. The solid line represents results of the fit and the 
dashed line shows the r/ c yield. 



C. B° — > A+p Branching Fraction Measurement 

From the B — > A+p fit yields in the B° —> ppK® and 
B + —> ppK* + modes given in Table IV we obtain the 
branching fractions shown in Table VIII. Averaging 
results for both modes and adding the errors in quadra- 
ture (except the systematic error on B— counting), 
we obtain the branching fraction B(B° — > A+p) = 
(21.0±^( S tat)±^(syst)±^(4+ B.F.)tH(B At ^ pK J)x 
1CP 6 . This measurement is consistent with the current 
value of B(B° -> A+p) = (21.9±|;| ± 3.2 ± 5.7) x lO" 6 
based on a single measurement by Belle [29] . 



B. B — > ( J/ip, r) c )h Branching Fraction 
Measurements 

Using the r\ c and J/ip yields from Table IV we obtain 
the branching fractions shown in Table VII. The values 
obtained are consistent with current world averages [12]. 
We also report the first evidence for the B + — » r\ c K* + 
decay, with a significance of 3.2<r (including systematic 
uncertainties). A sample of the maximum likelihood fit 
result for the charmonium region is shown in Fig. 2. 



D. B — > pph Charge Asymmetry Measurements 

The CP-violating charge asymmetry is defined as 
A ch = C% -N B )/(N W +N B ), where N B and Ng are 
the event yields in each of the categories of interest. The 
events yields are obtained from the maximum likelihood 
fit described in the Sec. IV, by integrating over the result- 
ing signal event weights for each of the two charge cate- 
gories separately. The resulting yields for all the modes 
(except B° — > ppK , which does not have information on 
the flavor of B meson) are summarized in Table IX. The 
measurements for the current modes are consistent with 
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TABLE VI: Summary of the experimental values for the branching fractions (xlO -6 ) of B — > pph and their comparison to 
two-body mesonic modes. The values for the two-body mesonic modes are taken from Ref. [12] unless otherwise noted. The 
values in bold are those presented in the current work. 



h 


Belle B -> pph [2] BARAR B -> pph 


B TT°h 


B -> p°h 


K+ 


6.0 ± 0.3 ± 0.4 6.7 ± 0.5 ± 0.4[1] 


12.1 ±0.8 




K° 


2.08±g;H ± 0.24 


3.0±0.5±0.3 


11.5 ±1.0 


5.4l?; 9 o [27] 


K*+ 


1n n+3.6 +1.3 
1U.J_2.8 -1.7 


5.3±1.5±1.3 


6.9 ± 2.4 


11.0 ±4.0 


K*° 


< 7.6, 90% CL 


1.5±0.5±0.4 < 3.5, 90% CL 5.6 ± 1.6 [28] 




1.68±g;i? ±0.12 


1.7±0.3±0.3 


5.5 ±0.6 


8.7 ± 1.1 



TABLE VII: Summary of the resulting branching fractions for the r\ c and J/tp modes (the order of the uncertainties is as follows: 
statistical, systematic, due to partial branching fraction correction where appropriate). The following values of branching 
fractions are used B{r) c -> pip) = (1.3 ± 0.4) x 10~ 3 and B(J/ip -> pp) = (2.12 ± 0.10) x 10~ 3 [12]. 



ppX 


Efficiency 


B(B^ Vc (pp)X)(W- 6 ) 


B(B -> Vc X) (10~ 3 ) 


Mode 


% 


Measured PDG [12] 


Measured PDG [12] 


ppK" 
ppK*° 
ppK*+ 


36.3 
23.7 
15.7 


0.83lo 2g ± 0.05 1.56±0.71 
1.03t° 2 2 l ±0-17 2.08±1.11 

i t-7+0. 56+0.45 
1 -° ' -0.46-0.36 


0.64+° 2o ± 0.04+°;^ 1.2±0.4 

0.80lg;^ ± 0.13±g;?| 1.6±0.7 
i 91+0.43+0.34+0.54 

1 - zl -0.35-0. 28-0.28 


Mode 


Efficiency B(B -> J/xjj{pp)X) (10~ 6 ) 


B(B -> J/tpX) (10 -3 ) 


ppK° 
ppK*° 
ppK*+ 


37.1 
25.0 
17.8 


1.871^6 ± °- 07 l-80±0.08 
2.82+°; 3 °+^ 2.78±0.20 
3.78lo;64lo.l 2.86±0.25 


0- 88 -o.i2 ± 0-°3 ± 0.04 0.85±0.04 
1.33+°;^ ±0.171^6 1.31±0.07 

1 78 +0.34+0. 13+0.09 1 or. n1n 
1 - 1 °-0. 30-0. 11-0.08 1-oOIEU.lU 



zero within less than three standard deviations. 



E. B(B° -f ©(1540)+p) Upper Limit Calculation 

As suggested in Ref. [3], we search for a pentaquark 
baryon candidate, 6 + , in the m pK o mass distribution of 

B° —> ppK° decays. If Q + decays strongly, there are 
only two possible decays modes: nK + and pK°. For this 
measurement we assume B(Q — > pKg) = 25%. From 
dedicated signal Monte Carlo we determine that the 9 + 
invariant mass resolution is represented by a sum of two 
Gaussian functions with common mean. The resolution 
of the main (secondary) Gaussian is 0.95(2.32) MeV/c 2 
and the wider Gaussian contributes 19% of the total. 
The overall resolution, defined as the full width at half 
maximum of the resolution function divided by 2.4, is 
lMeV/c 2 at the 6+ mass of 1.54GeV/c 2 . The 6+ pen- 
taquark signal efficiency is 30.8±0.1%. No events are ob- 
served in the + region of 1.52 < m pK a < 1.55GeV/c 2 
(see Ref. [30] for details). A Bayesian approach is used 
to calculate the upper limit of 0.20 x 10~ 6 at 90% con- 
fidence level, assuming Poisson-distributed events in the 
absence of background, including a multiplicative sys- 
tematic error of 7.1%. This value is consistent with and 
improves on the upper limit from the Belle Collaboration 
B(B° -> 9(1540)+p) < 0.92 x 10~ 6 [2]. 



F. Search for glueball /j(2220) in B — > pph decays 

We search for the narrow state /j(2220) by scanning 
through the 2.2 < m vp - < 2.4GeV/c 2 region with a 
30MeV/c 2 mass window in the final states ppK° and 
ppK* . This procedure is described in detail in Ref. [30]. 
The largest upper limits at 90% confidence level, includ- 
ing multiplicative systematic uncertainties of 2.7(4.3 and 
4.8)%, on the product of branching fractions are found 
to be B(B -> /j(2220)/i) x £(/j(2220) -> pp) < 4.5(7.7 
and 1.5) x 10~ 7 for h = K a (K*+ and K*°), respectively, 
assuming the /j(2220) width is less than 30MeV. 



VII. STUDY OF THE B -> pph DECAY 
DYNAMICS 

For decay dynamics studies, the maximum likeli- 
hood fit is performed using three variables (rriES, 
and Fisher discriminant) simultaneously over the whole 
Dalitz plot, with the exception of the ppTr + mode where 
we perform the fit in two regions: m pp - < 3.6 GeV/c 2 and 
m P p > 3.6 GeV/c 2 . 

The resulting background-subtracted efficiency- 
corrected Dalitz plots for all the modes are shown in 
Fig. 3. The main features of the Dalitz plots are expected 
to be the charmonium resonances (with J/ip — > pp and 
i] c — > pp bands most prominent), potential bands 
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TABLE VIII: Summary of the resulting branching fractions for B° —> A+p (the order of the uncertainties is as follows: 
statistical, systematic, the uncertainty in the ratio of the B(A~£ — > pK° / K*°) to the B(A£ — > pKn) [12] and the uncertainty in 
the B(A* — > pKn) value [12]). Note that efficiency does not include the sub-decay branching fractions. 



A+ decay Efficiency, % B(A+ -> pX)(l(T 3 )[12] B(B° -» yl+p)(l(T 6 ) statistical significance 
p~K° 25\4 23 ±6 15.1±|"g ± 3.2±1-|±|"1 3A 
pK*° 19.7 16±5 26.9t^ ! + 1 ^ _ + 3 4 1 °_ + 5 9 5 4 « 



TABLE IX: Summary of the asymmetry study: event yields 
for each of the categories of interest and the charge asymme- 
try. We observe no statistically significant CP asymmetries. 





Nb 


N B 


A ch 


Events Type 


S u 


ppK* y ' 


mode. 


Signal 


35±6 


28±5 


0.11±0.13±0.06 


Vc 


23±5 


13±4 


0.28±0.16±0.04 


JH> 


63±8 


43±7 


0.19±0.10±0.02 


Background 


15050±123 15443±125 


-0.013±0.006 




B+ 


— > ppK* + 


mode. 


Signal 


34±6 


18±4 


0.32±0.13±0.05 




7±3 


5±2 


0.20±0.28±0.07 


J/4> 


18±4 


16±4 


0.07±0.17±0.02 


Background 


5453±74 


5377±73 


0.007±0.010 




B+ 


— ► pp1T + 


mode. 


Signal 


97±10 


89±9 


0.04±0.07±0.04 


Background 


45434±213 44968±212 


0.005±0.003 



in B° — ► ppK° and B° — » ppK*° modes, as well as 
the low pp mass enhancements. Figure 4 shows the 
corresponding background Dalitz distributions. The 
combinatorial background events favor the edges of the 
Dalitz plot, as they are dominated by the inclusion of 
random low-momentum tracks. 

Background-subtracted efficiency-corrected m pp distri- 
butions are shown in Fig. 5 and are summarized in Ta- 
ble X. Although the m pp enhancement at low mass is 
quite prominent in the B° — ► ppK° and B + — > pp-n + 
modes, in the case of the B° — > ppK*° and B + — > ppK* + 
modes the statistics are too limited to draw a definite 
conclusion. Note that the shapes of the enhancement 
in B° — > ppK° and B + — > ppK + [1] are similar within 
the statistics of the measurements, in agreement with the 
theoretical predictions [5]. 

To shed light on the nature of this enhancement [10], its 
uniformity on the Dalitz plot has been tested. The Dalitz 
plot is divided along the diagonal m p h = m p h line and 
each of the two halves is projected onto the nearer axis. 
If the Dalitz plot is symmetric we expect the number 
of events in both projections to be the same. The corre- 
sponding background-subtracted efficiency-corrected dis- 
tributions for the signal events in all the modes are shown 
in Fig. 6. No asymmetry is expected to be introduced 
from variations in e mp _ which is charge-symmetric and 
slowly varying with m pp -. 

In the case of B° — > ppK°, there is no information on 
the flavor of the B meson and thus this study cannot be 
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FIG. 3: Dalitz plots for signal in the (a) 5° -> ppK°, (b) 
B+ -> pp-K + , (c) 5° -f ppK*° and (d) B + ppK*+ modes 
obtained using the weighting technique described in Ref. [25]. 
The positions of the following resonances are shown: r\ c in 
solid, J/tp in dashed and A c in dot-dashed lines. Note, that 
because of fluctuations and uncertainties, the signal rate in 
many bins is negative. The white areas correspond to regions 
with no entries. 



performed. For the B + — ► ppK* + mode there seems to 
be no difference between the two halves within the avail- 
able statistics. In the B° — ► ppK*° mode there might be 
a marginal excess at low m p h, which could be caused by 
the presence of the standard baryon resonances (such as 
A£), while in the B + — > ppn + mode there is a marginal 
excess of events at high m p h around 3.8 GeV/c 2 in the 
m p7T + half of the Dalitz plot, contrary to the result ob- 
served in the B + — ► ppK + mode [1]. No quantitative 
theoretical description of this correlation is available at 



present. Although the asymmetry in the low m pp band 
shown in [1] disfavors the possibility of the low mass pp 
enhancement originating only from the presence of a res- 
onance below threshold (such as the baryonium candi- 
date at 1835MeV/c 2 recently seen by BES [9]) the low 
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TABLE X: Summary of efficiency-corrected, background-subtracted number of events in m P p bins for B° — > ppK°, B + — > ppw 4 
B° — > ppK*°, and B + — > ppK* + modes. Note that contributions from the B decays to charmonium are removed. 



m P p range 


(GeV/c z ) 


B (i -> ppK K) 


B + — ► ppn + 


B" - ppK* u 


B + -> ppK* + 


1.876 - 


2.100 


160 ±50 ± 14 


96±30±9 


30 ± 37 ± 17 


280 ±109 ±39 


2.100 - 


2.350 


151 ±50 ±13 


74 ± 25 ± 9 


7±34± 16 


53 ± 90 ± 22 


2.350 - 


2.600 


80 ± 40 ± 7 


64 ± 24 ± 5 


3 ± 26 ± 5 


213±95 ± 18 


2.600 - 


2.850 


60 ±31 ±6 


83 ± 27 ± 8 


30 ± 26 ± 12 


189 ±86 ±22 


2.850 - 


3.150 


48 ± 23 ± 8 


22 ±21 ± 3 


32 ± 38 ± 5 


0±39±37 


3.150 - 


3.400 


-2 ±13 ±5 


37±18±3 


32 ±21 ±6 


112 ±66 ±33 


3.400 - 


3.600 


3 ±22 ±8 


-1±11±1 


64±27±8 


129 ±75 ± 14 


3.600 - 


3.900 


13±27± 5 


35 ±19 ±9 


47 ±31 ± 7 


84 ± 80 ± 17 


3.900 - 


4.150 


41 ± 25 ± 10 


11±9±1 


58 ± 30 ± 12 


30 ±65 ±9 


4.150 - 


4.387 


24 ± 25 ± 9 


-2±7± 1 


49 ± 30 ± 5 


105 ± 88 ± 17 


4.387 - 


4.580 


49 ± 29 ± 10 


-11±4±3 


-1±6± 1 


0±11±3 


4.580 - 


4.782 


25 ± 34 ± 6 


0±12± 17 


n/a 


n/a 


4.782 - 


5.139 


n/a 


12 ±18 ±19 


n/a 


n/a 
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FIG. 4: Dalitz plots for background in the (a) B° ppK° , 
(b) B + -> pp-K + , (c) B° -> ppK* {) and (d) B+ -> ppK*+ 
modes obtained using the weighting technique described in 
Ref. [25]. The expected positions of the following resonances 
are shown: rj c in solid, J/tp in dashed and A c in dot-dashed 
lines. The white areas correspond to regions with no entries. 
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FIG. 5: m P p distribution in data (points) for (a) B° — + ppK° , 
(b) B + -> pp-K + , (c) 5° -> ppK*° and (d) B + -» ppK* + 
modes obtained using the weighting technique described in 
Ref. [25]. The inner error bars show statistical uncertainty 
and the outer ones show statistical and systematic uncertain- 
ties added in quadrature. The histograms correspond to the 
relevant three-body phase-space signal Monte Carlo distribu- 
tions. The dashed line shows zero. The bin size is varying 
and is specified in Table X. Note that the contribution from 
the B decays to charmonium is removed. 



VIII. CONCLUSIONS 



statistics of the current modes does not allow to derive a 
definite conclusion. 



In summary, with 210 fb~ of data, we report evi- 
dence for the B° — > ppK*° decay with a branching frac- 
tion (1.5±0.5(stat)±0.4(syst))xlCT 6 , and provide im- 
proved measurements of branching fractions of the other 
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FIG. 6: m p h distribution obtained using the weighting tech- 
nique described in Ref. [25]: in red (open circles) for m v h > 
m-ph and in black (filled circles) for m p h < mph for (a) 
5° -> ppK°, (b) B+ -f PPTV+, (c) 5° -» ppK*° and (d) 
B + — > ppK* + . Solid line shows zero. Only statistical error 
bars are shown. 



£ — > pp/i modes, where h = ir + , K° and K* + . One key 
observation is that the pattern of decays for B —> pph dif- 
fers from that found for B — > tt°/i and B — > We also 
identify decays of the type £? — > X c -h — > pp/i, where /i = 
i^g, X* and if * + , and A C e = ?7c and J/ip. In particular, 
we report evidence for the B + — > rj c K* + decay with the 
branching fraction of B(B + — ► rj c K* + ) x S(^ c pp) = 
+0 56 (stat)l°;^(syst))xlO- 6 . We confirm the ear- 



lier observation of the _B° — > yl+p decay [29] and report 
measurements of the charge asymmetry consistent with 
zero in the B + — > pp7r+, B° — > ppK* a and — » ppK* + 
modes. No evidence is found for the pentaquark candi- 
date + in the mass range 1.52 to 1.55 GeV/c 2 , decaying 
into pK®, or the glueball candidate fj (2220) in the mass 
range 2.2 < m P p < 2.4GeV/c 2 , and branching fraction 
limits are established at the 10~ 7 level. 
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